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t^ ■ Abstract 

We present exact results for the periodic Anderson model for finite Hubbard 

^ . interaction < U < cxd on certain restricted domains of the model's phase 

o, 

, ^,1 diagram, in d = 1 dimension. Decomposing the Hamiltonian into positive 



semidefinite terms we find two quantum states to be ground state, an insu- 
lating and a metallic one. The ground state energy and several ground state 



\^ . expectation values are calculated. 

o 



I. INTRODUCTION 



'T3 I Since the last decades compounds with lanthanide or actinide elements as heavy-fermion 



systems have attracted great interest because of their interesting experimental behavior. 
Numerous strange phases are thought to be possible as the ground state of such a system, 

X" 

Vh ' ordered and disordered phases as well. As a theoretical background the periodic Anderson 

model (PAM) is one of the most important applicants to explain the physics in these systems. 
However, despite its relatively simple structure up to now the complete exact solution of the 
model is not known even in one dimension and there are mostly approximative results. In 
1993, using a method proposed by Brandt and GiesekusH, Stracka found the exact ground 
state in d = 1 for the strong coupling U = 00 case. His results were generalized for 
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d = 2 , 3 dimensions by Orlik and GulacsiS also in U = oo. This strong coupling limit is 
however an extremal case and is important to look for results in the physical domain where 
U < oo. Here we present for the first time exact results for PAM at finite U < oo, in 
d = 1 dimension. Two different solutions are reported, one insulating and a metallic one, 
both at 3/4 filling. The paper is organized as follows: The following Section is devoted to 
introduce the model and to show the positive decomposed form of the Hamiltonian used in 
the remaining part of the paper. Section III. contains the results of the calculations, and 
finally Section IV. concludes the results and closes the paper. 

II. THE HAMILTONIAN 

To start our study we use the Hamiltonian (H) of PAM as presented below: 

H = ff + t + Vo + V, + EfYl M,a + UU, (1) 

2,f7 

where T^ = J^i.a [^d c?j+i „ d'i,a + h.c] with d = c , f are hoppings for c and / electrons, 
Vo = Y.i,a [Vo ft,a Ci,<T + h.c] is the on-site and Vi = Y.i,a [ Vi ( cf^a fi+i,a + ft,a Q+i.a ) + h.c. ] 
is the first-neighbor hybridization term, U = J2i ^i ^ nf ^ is the Hubbard on-site interaction, 
and nf^ = df^ di^a ai'e the particle-number operators. 
To have a form of H appropriate for our goal let us transform Eq. (||) to the following form: 

^ = - E (i^p + bi' + kp + i^n + E «.^ <. + ( kp + bn E <. + 

i,cr i,a 

+ \v\' + Ef)Y.nl + f/E<T<i (2) 



«,(T 



where af^ = x cf^ + y c'l^i o- + ^ fi'a + "^ /i+i o- ^^^ x, y, z and v are given by the following 
equations: 

xy* = —tc, z V* = — tf , X z* + y V* = — Vq , xv* = zy* = —Vi. (3) 

In order to solve Eqs.(0) let us introduce m = y^. Using m, from Eqs.(0) we have z = mx, 
V = m* y and the square of x is determined by |xp = — -^^ ± y j^ ~ ^c t^- During 
the calculation we used the conditions 



tftc = V^, A + ^c>0, -^±\^,-tl—>Q. (4) 

^ 1 ' 4^*2 c - , 2m* V 4m*2 '^ m* " ^ ^ 

Going on with transforming H we get from Eq.(0) 

if = P + P' + (|x|2 + l^n ^ nl^ + {^\z\'' + \v\^ + Ef + U)Y. nl 

+ K - UL, (5) 

where L is the lattice size, P = J2i„ '^j.o- (^ta ^ P' = U J2i { ^ ^ ""-{t ~ ^,1 + ^{t '^U ) 
and K = — J2i,a ( kP + |yp ) ( |"^P + 1 )• After this step, from Eq.(^) in the case of 
Ef = (1 — |r7ip)(|sp + |yp) — U, the final form of H is given by 

H = P + P' + (|x|2 + \ynj2i<. + <J + i^ - f/L. (6) 

i,cr 

This form of H consists of constants K and UL, positive semidefinite operators P and P , 
and X^i.o-l^i'cr + ^{ct) which is related to the particle number as {n) = 1/L X]i,o-('^i.o- + ^{o-)- 
Because P and P are positive semidefinite, omitting them the ground state of the operator 
built up of the remaining terms must be a lower bound to the ground-state energy: 

El = {\x\^ + \y\^){n)L + K - UL. (7) 

An upper limit E^ for the ground state energy Eq can be derived from the variational 
principle, assuming some wave function | ■?/' ) to be the ground state. If there is any domain 
in the space of the model parameters with Ei = E^ = Eq then on that domain the exact 
ground state is | ^z^) and the exact ground state energy is Eq. This procedure will be followed 
in the next Section. 

III. THE OBTAINED SOLUTIONS 

In this Section we describe two different and exact ground-state solutions for PAM at 
finite U related to an insulating and a metallic phase. We start the presentation with the 
insulating case. 



A. The insulating state 



Consider that all the coupling constants are real. We have in this case from Eqs.( 
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where 02^ is an arbitrary phase. Inserting x and y into the expression of Ef given above we 



have E 

_ VqVi 
2tf 



f 



-il 



*/ A VoVi 



U. We also need the conditions tj tc = Vi , ^ > 4 



> 0, that are counterparts of those given in Eq.(H) but for the case of real coupling 



constants. In these conditions we have a lower limit for the ground state energy from Eq. (^ 
and Eq.(|[), in the case of ( n ) = 3: 
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tf - vr 

To get also an upper bound we call the variational principle. Let the trial wave function be 

iV'i) = n(<T + "^/^)(<i + m/+)(«,c+ + Ac+^ + 7./+T + S^Kl)\0) (10) 

i 
i 

where e^ and z/j are related to the quantities ai , Pi , ji , 6i. Since | ipi ) sets three electrons on 
any lattice site it belongs to band filling (n) = 3. Using Eq.(|lOD we have an upper limit for 
the ground state energy: Eu,i = {ipi\H\tpi) / {ipi\ipi). Since P | ■j/'i ) = P \ipi) = 0, it follows 
that Eu,i = Ei^i = Eq, where Eq is the ground state energy. Collecting all the conditions 
weused: f/> 0,t;t, = ^1^^ > 4, -^ > 0, E^ = - ( 1 - i^) 1^ - f/, (n) = 3, 
the ground state is given by Eq.(|TOD and the ground state energy is of the form of Eq.(^. 
To get a deeper insight into the nature of the state represented by \ ipi) w^ show some 



expectation values calculated with it. Since | ipi ) sets three electrons on each site it is 
obvious that any way of moving an electron from a certain site i to an other j ( 7^ i) results 
an orthogonal state to it so we have 

(T,) = {Tf) = {V,) = 0. (11) 

Expectation values of other terms after normalization are given as 

, f . m^ -^ 

(<t) = TTz::^. + 



■T / i + rn^ ^ {ef + iyf){l + m^)' 

( n-' ) ^ -I- - 

^ ''-^' 1 + m2 ^ ( ef + z/2 ) ( 1 + m2 ) ' 

\ ^i,T ) - YT^ ^ ( ef + z/2 ) ( 1 + m2 ) ' 
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It is known that in the case of large U the periodic Anderson model goes to the Kondo lattice 
model which has an interesting insulating ground state at half-filling called the Kondo- 
insulator. As can be seen, PAM in strong coupling limit is able to give an insulating phase 
not only at half-filling, but also at 3/4 filling as presented here. From Eg . (|I^) we see that 
I ipi ) contains all the possible spin directions (i.e. is degenerate in spin) and thus describes 
an insulating paramagnetic state with zero net magnetic moment. Figure 1. and 2. present 
the emergence domains of | ^1 ) in the parameter space. 

B. The metallic state 

Let consider the parameters tc and tf real, Vq and Vi imaginary, the Hamiltonian re- 
maining of course Hermitian. We have from Eq. (^ 
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where (^^ is arbitrary. Similarly to the previous case, from x and y given here we have 
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ypM/ 4*2 I tc~ U- The only condition coming from Eq.(|) is tjtc = V^ 
because the others are automatically valid. In the case oi (n) =3 from Eq.(|^) and Eq.(|12|) 
we find a lower limit for the ground state energy as 
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The upper limit to Eq is derived with the trial function 

where «« , /3j arbitrary, x/"^ = df^ or e/^^. V i and 
<<x = a: ( c+^ + m /+^ ) , e+^ = y ( c+ ^^^ + m /+ ^^^ ) 

and X, y are given in Eqs.(|T^. It can be verified that | ip2 ) belongs to (n) = 3. The upper 
limit to the ground state energy is £'m,2 = (V'2|-f^|'^2)/(V'2 |'^2)- Since P 1 7/^2 ) = P |V^2) = 0, 
it follows that £'^,2 = -E«,2 = Eq is the ground state energy. Collecting all conditions: U > 0, 

j2 /t/2i/2 

/ = —2(1 — j^ ) J -jjT~ +t1 — U, {n) = 3 we find that the ground 



tj tc 
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state is given in Eq.(|T^) and the ground state energy is as shown in Eq.(|T3|). In order to 
understand the physical behavior of this state let us examine expectation values calculated 
with it. Imposing the constraint in Eq . (pjlf) a = Oi, j3 = /3j Vz we can express | ?/'2 ) in the 
momentum space: 



^-2 ) = n [ ^ ( 4t + ^ /i!"T ) + ^ ( 4+i,T + "^*//+i,T ) ] X 



X [^ (4i + ^fu) + y(4+u + ^* /iVu ) ] [« /i!^ + /5/i^] 1 0) 
= n [ E <,T ^''''' ] [ E 4„ie"^'^^] [ E ( « /i,T + /^ /.ti ) ^''''' ] 1 ) 

j fcl fc2 fc3 

with 4_, = ( X + 1/ e-''^' ) c+^, + ( a: m + y m* e'^'^^ ) /+, , (15) 

From this expression we find 

k 

= M Y[[ela cl^ cl^ f^^ + d 62 a 4| /fc+| f^^ + 

+ e? /5 c+^ c+^ /+^ + ei 62 /5 /,+T c+^ /+^ ] I 0) 

ei = X + y e~''' , el = \x\^ + \y\^ + 2tcCos{k), 

62 = X m + ym* e~^^ , 63 = |m|^ ( k^ + ll/T ) + 2 ^^ tc cos{k) , 

JVf = V (-1)^ g-3i[lp(l)+2p{2)+...+iVp(7V)] /-^gN 

{^) 
where P denotes a permutation (p(l) , p(2) . . ■p{N) ) of ( 1 ... A^ ). 
The norm of \1IJ2) is given as 

(^21^2) = \K\' ( \a\' + \f3\' ) n ( |ei|' + \ei\' \e2\' ) 

k 

= |K|2(|a|2 + l/^n n["^(l + l^n + ^tcUcos{k) 

k 

+ 4tc(l - |mp)cos^(A;)] , 

where m = |xp + ||/p. After this step normed expectation values of different quantities can 
be given: 

|/5p u^ + 4tcMC0s(A;) + 4t2cos2(A;) 



^"-^(^)^ = lap + I/3P 



a|2 + |/3|2 m2 ( 1 + |^|2 ) + 4tcU cos{k) + 4 tc ( 1 - I^P ) cos^k) ' 



+ 



ap M^ + 4tcMcos(A;) + 4t^cos2(A;) 



|a|2 + |/5|2 m2(i + |^|2) + 4tc'«cos(A;) + 4 t^ ( 1 - |m|2 ) cos2(A;) ^ 



lap 
^"^'^^"^^^ = kP + I^P 
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|/3|2 |m|2(M2 - Atlcos^ik)) 



<-Ai(^)) = |«p + l.p 



aP + |/?|2 u'^{l + |m|2 ) + AtcUcos{k) + 4 tc ( 1 - |mp ) cos2(A;) ' 



|ap |mp (m^ - Atlcos'^{k) ) 
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|a|2 + |/5|2 ^2(1 + |m|2) + At^ucos{k) + 4 te ( 1 - |m|2 ) cos2(A;) 

-|/3|2 ( u + 2 te cos(fc) ) ( K) + 2i Vi sm(A;)) 
^2 ( 1 + |m|2 ) + 4 tc w cos(/c) + 4 tc ( 1 — I^^P ) cos2(/c) ' 

-\a\^{u + 2 tc cos(A;) ) {Vo + 2iVi sin{k) ) 
^2(1 + |m|2) + 4tcUcos{k) + 4tc(l - |m|2 ) cos2(fc) 



^^^■. 



fc 
From Eqs.fp!?!) it is seen that | 'i/'2 ) describes an insulating state. This state is also para- 
magnetic because | '?/'2 ) is degenerate in spin. This degeneration is given by the arbitrary 
quantities a and /3. Furthermore, because (rik) = {ric^i^k) + rac,j,(/c) +n/^f(/c) +^/,j,(fc)) = 3 
is k independent, in (uk) there is no any change present at the Fermi-momentum kp- Fig 
3. presents the | iIj2 ) phase in the parameter space while Fig 4. shows the relative situation 
of the phases | ipi ) and | V^2 ) • 

C. Comparison of the two states 

Let us compare the wave functions | ipi ) and \ 1^2)- The mathematical relationship 
between Eq. (|T0|) and Eq. (0) can be seen in more transparent form after an orthogonalisation 
procedure. For this reason we decompose Eq.(p^ as follows 

I 4) = \^2) - 2n<T<i(«^/^T + A/i)|0) 

i 

= \^2) - 2n<T<i(«^/^ + A/+i)|0) (18) 

i 

The wave- vector | ^2 ) is zero if m is a real number. This is seen from the fact that 1^2) 
contains at least one term of the form 

ej> d^+i,a = xy{ C++1 ,, + m* //^^^^ ) ( c++i ,, + m /^+ ^ ,, ) (19) 

which, due to the fermionic anticommutation rules, is zero if m is real. It is also seen that in 
the case of infinite lattice Jli c^j^T ^ti ( "« /vr + ^^ ^U ) I "^ ) = H j e^-^ e+| ( a* /i^| + A /j|| ) 1 ) = \ipi). 
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The state | ipi ) sets rigorously 3 electrons on each site, while | ip2 ) contains at least one site 
with 4 electron an one with 2 electrons. As a consequence, | ipi ) and | ■i/'g ) ^^^ orthogonal. If 
m is a real number then the ground-state is {ipi), and, if m is imaginary then the ground- 
state is I '?/'2 ) = I ^"2 ) + 2 I V^i ) . The fact that {ipilTclipi) = but {4'2\Tc\ip2) ¥" is 
attributed to the effect of | '?/'2 ) i^ the ground state. 

IV. CONCLUSION 

Using a positive semidefinite decomposition of the Hamiltonian we have found exact 
results for the periodic Anderson model at finite U in d = 1 dimension and 3/4 band filling. 
The result detailed in first is for the case when all the coupling parameters have real value. 
In this case we found that the ground state is Eq. ([To|) and the ground state energy is of 
the form of Eq.(||). The ground state Eq. ([lOD describes a paramagnetic insulating state 
that interestingly is present also in the strong coupling limit (i.e. Kondo lattice case) at 
three-quarters band filling. 

The second solution presented is for the case when the hopping amplitudes are real and 
the hybridization amplitudes are imaginary (the Hamiltonian remaining Hermitian). We 
found a ground state presented in Eq.([l^) and a ground state energy as given in Eq.(|l^). 
This state has nonzero expectation values for the c and / hoppings so it is a metallic state. 
We have found that the state has no discontinuities in n^ at kp. 
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FIGURE CAPTIONS 






Fig.l.: The I ipi ) insulating state in the phase diagram. Notations are X = j^,Y 
Z = J-, a = -j^. Setting for example a = — 0.5 the quantum state \ipi) \s ground state 
when Z = — ( 1 — ^) XY — a^ with the constraints— X F > 0, ^ > 4. 

Fig. 2.: A part of the | ipi ) phase surface. Notations are the same as in Fig.l. 

Fig. 3.: The | ■?/'2 ) metallic state in the phase diagram. For this Figure the notations are 

T/2 '1/2 jj rp 

X = -T, Y = -§-, Z = j^, a = 7^. Goosing a = — 0.5 quantum state I ■?/'2 ) is found to 
*/ '/ '-/ '-/ 

be ground state when Z = — {1 — -^) JX^ + ^XY — a. 

Fig. 4.: The phases | ipi ) and | ■?/'2 ) in the phase diagram presented together. 
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FIGURES 




